Submitted to ApJ, October 10, 2005 

Preprint typeset using I4T]^X style emulateapj v. 6/22/04 



HOW MUCH MASS DO SUPERMASSIVE BLACK HOLES EAT IN THEIR OLD AGE? 
Philip F. Hopkins 1 , Ramesh Narayan 1 , & Lars Hernquist 1 

Submitted to ApJ, October 10, 2005 

ABSTRACT 

We consider the distribution of local supermassive black hole Eddington ratios and accretion rates, account- 
ing for the dependence of radiative efficiency and bolometric corrections on the accretion rate. We find that 
black hole mass growth, both of the integrated mass density and the masses of most individual objects, must 
be dominated by an earlier, radiatively efficient, high accretion rate stage, and not by the radiatively inefficient 
low accretion rate phase in which most local supermassive black holes are currently observed. This conclusion 
is particularly true of supermassive black holes in elliptical host galaxies, as expected if they have undergone 
merger activity in the past which would fuel quasar activity and rapid growth. We discuss models of the time 
evolution of accretion rates and show that they all predict significant mass growth in a prior radiatively efficient 
state. The only way to avoid this conclusion is through careful fine-tuning of the accretion/quasar timescale to 
a value that is inconsistent with observations. Our results agree with a wide range of observational inferences 
drawn from the quasar luminosity function and X-ray background synthesis models, but our approach has the 
virtue of being independent of the modeling of source populations. Models in which black holes spend the 
great majority of their time in low accretion rate phases are thus completely consistent both with observations 
implying mass gain in relatively short, high accretion rate phases and with the local distribution of accretion 
rates. 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 



1. INTRODUCTION 

The growth of quasars and supermassive black holes is a 
topic of fundamental interest in cosmology, critically inform- 
ing theories of black hole formation and accretion models, and 
the role of black hole energetics in the X-ray, UV, and infrared 
backgrounds and their consequences for reionization. It is 
now believed that supermassive blac k holes reside at the cen- 
ters of most if not all galaxies (e.g. iKo rmendv & Richstone 
119951 iRichstone et 10119981 fKormendv & Gebh ardtl l200lT 
Furthermore, recent discoveries of correlations between 
masses of black holes in nearby galaxies and either the mass 
(Magorrian et al. 1998; McLure & Dunlop 2002; Marconi 
& Hunt 2003) or velocity dispersion (i.e. the M-a relation: 
Ferrarese & Merritt 2000; Gebhardt et al. 2000) of spheroids 
demonstrate a fundamental link between the growth of super- 
massive black holes and galaxy formation. 

Observations of the nearby Universe suggest that rapid 
black hole growth may be related to gas inflows in 
the center s of galaxies, powe ring intense infrared emis- 
sion (e.g. Sanders et al. 1988), with mergers or interac- 
tions of galaxies providing the gravitational torques nec- 
essary to drive large quantities of gas to galaxy centers 
(e.g. Barnes & Hernquist 199L 1996) to fuel both star- 
bursts llMihos & Hernquistl 119941 Il996|) and rapid , self- 
regulating black hole growth (Di Matt eo et al.l 120051) . This 
picture of quasar activity, whether represente d in a nalyti- 
cal or numerical models of the M-a relat ion (Silk & Rees 
1998; Fabianll 19991 IcioTt i & Ostrike^ l200lt IDi Matteo etafl 
2005), semi-analytical models or computer simulations 
of galaxy and quasar forma tion (IKa uffmann & Haehnelt 
2000; iHaiman & Menoull2000t IDi Matteo et al.ll2003t 120041 
Wvith e & Loebl 120031 iGranato et al] 120041) . or modeling of 
quasar light curves and the quasar and elliptical galaxy lumi- 
nosity functions (Hopkins et al. 2005a-g), implies that black 
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hole growth should be dominated by radiatively efficient, rel- 
atively short phases at high accretion rates (> 10% of Edding- 
ton). 

This picture is supported by inferred black hole accre- 
tion histories from the qu asar lu minosity function (e.g . 
I Soltanl 119821 ISalucci et all I1999L lYu & Tremaiiiel [ 20021 
Marconi et al] 120041 IShankar et all 120041 iMerlonil boo4 
Merloni et al. 2004) and synthesis mod els of the X-ray back- 
ground (e.g. Comast ri et alll995llGilli et alJl99^|Elvis et alJ 
l200l lUeda et alJ 120031 ICaol l2005t iBarger et alJ 120051) . as 
well as direct obs ervations of qua s ar accretion rates at 
diffe rent redshifts dVjgtergaard 20041 iHeckman et alJ 12004 
iMcLure & Dunlopl 120041) or models which combine all of 
these based on quasar evolution an d the observed luminos- 
ity function (Hopki ns et al J 12005 del) . However, these infer- 
ences are generally model dependent, and more importantly 
are based on quasar observations which, while in some cases 
probing X-ray luminosities well below L* and much fainter 
than bright optically selected quasars, do not include the (at 
least locally) very large low accretion rate, radiatively inef- 
ficient population. It is not clear how strongly constrained 
is the possibility of black holes gaining a significant amount 
of mass in much lower accretion rate phases after periods of 
early, bright activity; in these latter phases the radiative effi- 
ciency may drop, and thus the black holes will not usually be 
observed. 

Various authors have attempted to use observational esti- 
mates of the distribution of quasar accretion rates as a func- 
tion of redshift to answer the question of whether or not 
all black holes pass through and gain their mass in a high- 
Eddin gton ratio bright quasar phase (e.g.. IHeckman et ail 
120041 IMcLure & Dunlopll2004t iKollmeier et alJl2005l) . Al- 
though the determination of the distribution and evolution of 
quasar accretion rates remains a key outstanding problem, se- 
lection effects limit the questions these high-redshift observa- 
tions can answer. These surveys are generally complete only 
down to accretion rates ^0.1 of Eddington, and at any given 
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redshift only a small fraction of black holes are "active" with 
Eddington ratios above this fraction. Surveys which can probe 
lower accretion rates have yielded very different results for 
the number and distribution of objects with accretion rates in 
the range ~ 0.01-0.1 Eddington (compare, e.g. the estimates 
of Vanden Berk et al. 2005 from SDSS and Kollmeier et al. 
2005 from AGES). Consequently, although the observations 
described above may suggest otherwise, these high-redshift 
measurements of the accretion rate distribution at large Ed- 
dington ratios cannot determine whether or not a significant 
population of black holes exists at accretion rates ~ 1% of 
Eddington, sufficient to dominate their mass growth and total 
mass density by z = 0. 

Furthermore, estimating the mass growth of individual 
black holes (as opposed to the integrated growth of the popu- 
lation as a whole) in the observed high-accretion rate phases 
is not directly possible, as the duration of this phase is degen- 
erate with the rate at which quasars are "triggered" or "ac- 
tivated" at these rates. Even indirectly, neither quantity is 
well-constrained, as e.g. the duration of this phase is currently 
restricted obse rvationally only to the range ~ 10 6 - 10 8 yr 
l lMartinil2004l) . 

At low accretion rates (relative to Eddington), black 
holes transition to radiatively inefficient accretion flows 
(RIAF/ADAF e.g.. iNaravan & Yil Il995bt lOuataertl EOQlt 
Naravan 2004; Yuan & Naravan 2004), making these objects 
especially difficult to observe, and implying that observational 
surveys may be an order of magnitude less deep in terms of 
intrinsic accretion rate than they are in terms of a radiative 
(L/LEdd) Eddington ratio. It is therefore especially important 
to consider the distribution of accretion rates at low redshifts, 
where objects with low Eddington ratios can be reliably de- 
tected, as a constraint on possible black hole growth in low 
Eddington ratio states. The fact that the population of low 
Eddington ratio sources, which must exist at all moderate and 
low redshifts (owing to limits on the integrated black hole 
mass and X-ray luminosity density as discussed above), is ob- 
served at z = to extend to accretion rates as low as 10~ 6 of 
Eddington strongly emphasizes this point, as these accretion 
rates are unlikely to be observable in high-redshift surveys in 
the near future. 

In this paper, we consider the z = observed distribution of 
accretion rates and determine whether or not it allows for sig- 
nificant mass growth in radiatively inefficient, low accretion 
rate states. There are three basic questions one could ask: 

(1) In which state — radiatively inefficient, low accretion 
rate, or radiatively efficient, high accretion rate — do black 
holes spend most of their time! 

(2) In which state do black holes gain most of their massl 

(3) In which state do black holes radiate most of their lumi- 
nosity! 

According to current wisdom, the answer to question 1 is 
the radiatively inefficient state, and that to question 3 is the 
radiatively efficient state. However, question 2 is rather del- 
icate; it is by no means straightforward to deduce its answer 
from the other answers. The purpose of the present paper is 
to investigate question 2 in depth using the z = distribution 
of black hole masses and accretion rates. 

We describe th e observed sample from | Hol (l2002) and com- 
pare with that of iMarchesini et al. (2004! i, and discuss possi- 
ble selection effects in § 12.11 In § 12. 21 we discuss the depen- 
dence of bolometric corrections and radiative efficiency on 
accretion rate and our accounting f or th is in determining the 
distribution of accretion rates. In § I2.3l we consider the dis- 



tribution of accretion rates as a function of host galaxy mor- 
phology. We discuss several models for the mass growth of 
black holes and their implications for the possibility of mass 
gain in low accretion rate states in §|3] and finally summarize 
our conclusions in §0] 

2. THE OBSERVED DISTRIBUTION OF LOW-EFFICIENCY 
ACCRETION RATES 

2.1. The Sample and Possible Selection Effects 

We consider the observed sample of Ho ( 2002), itse lf an up - 
date of the compilation of black hole masses in iHol (11999a). 
which includes black hole masses, B-band luminosities, red- 
shifts, and host galaxy morphological information for 80 
nearby supermassive black holes. The measurement methods 
for black hole masses and luminosities, and possible selection 
effects, are discussed in detail in IHol (|2002), but we briefly 
review them here. 

The majority of the black hole masses are derived from spa- 
tially resolved observations of gas and/or stellar kinematics 
(for a review of the methods, see e.g. Kormendy & Rich- 
stone 1995). There are a few exceptions, namely Sgr A* in 
the Milky Way, whose mass is estimated based on prop er mo- 
tions of individual stars ("e.g JSchodel et alT20 02: Ecka rt etail 
20021 Idhez et al.l2 005 ), four galaxies which have strong wa- 
ter maser emission (NGC 1068, NGC 4945, Circinus, NGC 
4258; Greenhill et al. 1996, 1997, 2000; Miyoshi et al. 1995, 
respectively), and Arp 102B for which the orbital pe riod o f 
th e accretion disk has been determined (tNewman et alJl997h . 
In IHol ( 120021) . the differences between masses estimated by 
different means of dynamical modeling of the integrated spec- 
troscopy (primarily HST spectroscopy of stellar dynamics in 
the central regions of each galaxy) are discussed. Three- 
integral dynamical modeling of the galaxy surface brightness 
profiles and stellar absorption lines is used to estimate the 
black hole masses, which tends to give slightly lower black 
hole masses th an the simplified axisymmetric two-in tegral 
modeling of e.g. Magorrian et al. ( 1998) (van der MarelH 999: 
Ferrarese & Merritt 2000; Gebhar dt et al.l2000l) . 

For more luminous or distant active galactic nuclei (AGN), 
direct dynamical measurements are no longer feasible, as the 
AGN luminosity will overpower the central stellar emission 
features. Thus, a large fraction of th e sample masses are deter- 
mined via reverberation mapping (iBlandford & McKed 19821) 
of the central emi ssion lines. T h e un certaintie s of this m ethod 
are dis cussed in | Ho | (Il999al): IWa ndel. Peterson. & Malkan 
( 1999): lKaspi et all (120001) : lMcLure & Dunlop|(l2001l) : lKrolikl 
(2001), but it is generally true that black hole masses es- 
timated via reverberation mapping agree with those deter- 
mined from dynamical estimates to within a fac tor of ~ 2- 3. 
The r everberation-mapp ed o bject s oflHo ( 2002) are compiled 
from lKaspi etalJ d2000h and lHol lll999al) . in which the indi- 
vidual objects are described further. 

The optical (B-band) luminosities of the sources in the sam- 
ple are determined using the correlation between H/3 lumi- 
nosity and B-band abso lute magnitude dYeel[l980t [Shuder 
119811) . as calibrated by Ho & PengJ ( 120011) . which is in prin- 
ciple a more isotropic quantity than optical continuum lumi- 
nosity and is measurable in obscured (Type II) AGN. Where 
the B-band magnitude is directly available (including all of 
the reverberation-mapped sources), it agrees reasonably well 
with this determin ation, wit h no significant systematic bias 
llHo & Pengll2001t|Holl2002l) . We do not consider the 17 PG 
quasars in the Ho (2002) sample, as these objects are chosen 
from yet another parent sample, at non-negligible redshifts 
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FIG. 1. — The distribution of bolometric luminosities (left) and radiative 
Eddington ratios (A = L/Lrjdd. right) in the Ho 12002) sample (solid his- 
tograms) and the Marchesini et al. 1 2004 ) radio-selected sample (dot-dashed 
histograms). Open histogram shows the cumulative sample of He 1 2002 ) , 
horizontally shaded shows the distribution for objects with stellar or gas dy- 
namical mass determinations, vertically shaded shows the distribution for ob- 
jects with reverberation mapping mass determinations. 



z = 0. 1 -0.2, introducing stronger selection effects. Mo reover, 
they a re all at reasonably large Eddington ratios dHolll999al 
2002), and thus are members of the class of higher-redshift, 
high-accretion rate radiatively efficient objects in which we 
are explicitly not interested. We have however checked that 
including them in our total sample results in a negligible dif- 
ference in the cumulative Eddington ratio distribution and our 
subsequent calculations, after normalizing the relative frac- 
tion for the effective volume of the subsample (i.e. converting 
the observed number to a number density) in Figure[2]below. 

We briefly consider possible selection effects of the differ- 
ent methods by which the black hole masses in this sam- 
ple are measured. As noted above, strong or bright AGN 
will overwhelm the stellar emission in the central regions of 
the observed galaxy and make direct measurements of stel- 
lar kinematics difficult. Therefore, the objects with masses 
measured via stellar kinematics may be systematically bi- 
ased towards lower nuclear luminosities, and correspondingly 
lower Eddington ratios and accretion rates. However, the 
reverberation-mapped objects are expected to have the oppo- 
site bias, as this method of mass measurement requires resolv- 
ing the AGN emission lines, and thus may bias the observed 
sample to higher Eddington ratios and accretion rates. 

Figure [2 shows these effects explicitly. Here, the distri- 
bution of bolometric luminosities of our entire adopted Ho 
(2002) sample is shown (open histograms, left panel). Since 
we simply wish to show the luminosity range sampled by 
the observations we use a constant bolometric correction (L = 
1 1 Lb ), but we adopt more accurate bolometric corrections in 
§ 12.21 The right panel shows the corresponding distribution of 
radiative Eddington ratios. Horizontally shaded histograms 
show the distribution for objects with masses determined us- 
ing gas or stellar dynamics, whereas vertically shaded his- 
tograms show the distribution for objects with masses deter- 
mined using reverberation mapping. There is a considerable 
gap in the luminosities probed around 10 9 L©, and a simi- 
lar gap in Eddington-scaled luminosities around 10~ 4 - 10~ 2 . 
Owing to the minimal overlap in the bolometric luminosities 
covered by either method, it is difficult to reliably determine 



whether the observed gap in the Eddington ratio distribution 
is real or represents the combination of the selection effects of 
the two methods. 

To test this, we consider the sam ple of black hole masse s 
and bolometric luminosities from Marches ini et alJ (2004), 
shown in both panels as dot-dashed histograms. We consider 
the FR I sample studied by these authors, based on the 298 
radio galaxies and 53 radio loud quasars of the 3CR catalogue 
( Spinrad et al. 1985). The black hole masses are estimated 
from the stellar luminosities of the host galaxy bulges (ex- 
cept for the bright radio loud quasars, which we do not con- 
sider here for the same reason we neglect the high accretion 
rate quasars in the Ho [2002] sample, although again we have 
checked and find that the relative number density is not large 
enough to change our results). The bolometric luminosities 
are based on fitting the nuclear optical luminosities to the light 
profiles of the observed galaxies. 

The advantage of considering the Marchesini et al. (2004) 
sample is clear in the left panel of Figure [H as it cove rs the 
range of bolometric luminosities in which the Ho (2002J) sam- 
ples show a deficit, yielding a combined sample with a smooth 
(non-bimodal) luminosity distribution. The Eddington ratio 
distribution for these objects is shown in the right panel, and 
although it peaks at a higher Eddington ratio (~ 1 -4 ) t han 
the stellar/gas dynamics-measured sample of iHol J2002I) . it 
has a similar qualitative behavior. Furthermore, this higher 
peak most likely owes to the higher luminosity limit in the 
iMarchesini et alJ J2004I) sample. 

Although including the lMarchesini et all (|2004) sample fills 
the gap in bolometric luminosity, the gap in Eddington ra- 
tio at 10~ 3 - 10~ 2 is pre served, a point which is discussed in 
Marc hesini et alJ (|2004l). I f we consider the combined distri- 
bution fromthelHgy^^) an d Marche sini et alJ J2004I) sam- 
ples, the low-Eddington ratio distribution peaks at roughly 
10~ 5 , more than low enough that our subsequent results are 
all qualitatively unchanged regard less of whether o r not we 
include the objects in the Marchesini et al. (2004) sample. 
Based on this, the "gap" at Eddington ratios ~ 10~ 2 does not 
appear to be a selection effect, and we can rely on the Ho 
(2002) sample in our following analysis. 

In the following, we ge nerally do not include the sample 
of Marchesini et al. (2004) except where otherwise specified, 
because it is unclear how to correct for possible bias intro- 
duced in considering a radio-loud sample, given the correla- 
tion betwe en Eddington ratio and radio loudness observed in 
IHol J2002I) . For further discussion of these distinctions be- 
tween samples and possible sel e ction effects, we refer to lHol 
J2002I) : IMarchesini et all j2004l) : ljeste3 120051) . 

2.2. Determining the Accretion Rate 

The bolometric luminosity, which we denote by L, of a 
black hole accreting at a rate M is L = e r Mc 2 , where e r is the 
radiative efficiency. Therefore, the Eddington rate of accre- 
tion for a black hole of mass M is M E dd = L^iM) / e r c 2 , where 
the Eddington luminosity of a black hole is well-defined, 

L Edd « 3.3 x 1O 4 L (M/M Q ). (1) 

Because the Eddington accretion rate depends on the radiative 
efficiency e r , which itself can depend on the accretion rate, 
these definitions are circular unless we define the Eddington 
accretion rate with respect to a canonical radiative efficiency 
e°. We follow standard practice (e.g. Narayan & Yi 1995b; 
Esin, McClintock, & Narayan 1997; Jester 2005, but see also 
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Chen et al. 1995; Blandford & Begelman 1999; Marchesini et 
al. 2004 who choose e° r = 1) in defining 



e? = 0.1, 



(2) 



which corresponds to the typical radiative efficiency 
expected from accretion through a standard efficient 
iShakura & Sunvaevl (119731) thin disk. We then define the Ed- 
dington rate as 



M Edd : 



^Edd(M) _ M 

e°c 2 t s ' 



where 



t s = A2 x 10 yr 



(3) 



(4) 



is the Salpeter time, the mass e-folding time for a black hol e 
accreting at the Eddington rate with e r = 0.1 (Salpeter 1964). 
We can further define the dimensionless luminosity Eddington 
ratio A and mass accretion rate Eddington ratio rh by 



A=L/L Edd , m=M/M Edd . 
With these definitions, we have 



A = 



e° r Mc 2 /t s 



. e r (m) 
= m — t[ — , 



(5) 



(6) 



where we allow explicitly for the possibility that the radiative 
efficiency may depend on rh. 

We can also define the black hole growth timescale, which 
corresponds to the mass e-folding time for growth at a con- 
stant Eddington ratio, 



M 1 

fBH = 7v = — ts- 

M m 



(7) 



Note that this further lets us define a special value of rh, 

m H EEf s /fHubble ~ 0.003, (8) 

which is the dimensionless accretion rate for which the growth 
time is equal to the Hubble time (~44Gyr). For significantly 
lower accretion rates rh <§C /«h, the growth timescale is much 
greater than the age of the Universe, and a given black hole 
cannot possibly gain a significant fraction of its mass at that 
rh. 

For a simp le estimate of the accretion luminosity L, we 
have followed Ho (2002) and converted from B-band to bolo- 
metric luminosities adopting a constant bolometric conver- 
sion (L = cbLb, cb = 17). However, the actual bolometric 
conversion varies, rangi ng from ~ 11-17 for bright AGN 
and quasars (e .g. lElvis et alj fl994). to - 24 (e.g. lHolfT999bT: 
iHo et alJ2"0 00) for the most low-luminosity RIAF/ADAF sys- 
tems. For bright AGN and quasars, the bolometric conver- 
sion cb appears to depend just on luminosity, although this 
of course may instead reflect a dependence on the accre- 
tion rate. For objects with moderate to large accretion rates, 
we adopt the lu minosity-dependent bolometric conversions of 
iMarconi et alJ |2004), based on observations of the quasar 
spectrum over a wide range of wa velength s and as a func- 
tion of luminosity ("e.g.. lElvis et al.| |l994|:|George et ail|1998|; 
I Vanden Berk et alJ200UlPerola et alJ2002UTelfer et alJhooi 
llJeda et alJ2003HVignali et al.l2003h . This gives a B-band lu- 
minosity 



log(L B ) = 0.80-0.067£ + 0.017£ 2 -0.0023£ 3 



(9) 



where C = log(L/L Q )- 12. The important point is that c# 
decreases with increasing luminosity; i.e. the brightest ob- 
jects are the most dominated by the optical-UV portion of the 



spectrum (see also, e.g. Wilke s"et al.ll994tlGreen et all 1993 
Stratevaetall2005l). 

At low luminosities, we consider th e sub- sample of low- 
luminosity AGN from IHoI dl999bl) and lHo et al] (I2000I) . for 
which bolometric luminosities from radio, optical, UV, and 
X-ray observations have been determined directly. Calculat- 
ing the cb correction for each of these objects, we find no evi- 
dence for a direct correlation with L, but we do find marginal 
evidence for a correlation of eg with observed Eddington ra- 
tio A = L/L Edd . We can fit this dependence roughly to a power 
law and obtain the best fit 



c B w 11 (A/0.01) 



-0.25 



(10) 



As essentially all these objects have rh < rh cr n (defined below, 
but roughly m cr it ~ 0.01-0.1), we do not need to worry about 
a "break" in the slope in terms of A induced by the change in 
dependence of e r on A (i.e. on m) at m CI i t . Again, we find that 
the ratio of total to B-band luminosity decreases with increas- 
ing luminosity (though in this case with dimensionless lumi- 
nosity). This is also consistent with theor etical spectral mod - 
els of advection-dominated accretion (e.g. Mahadevan 1997), 
which generally feature an increasingly hard X-ray spectrum 
and possibly increased contribution from jets at low rh. 

In order to determine which objects should receive which 
bolometric correction, we divide our sample in two on the 
basis of our original calculation of A. For objects with 
A > m cs it (defin ed below), we apply the bolometric correc- 
tions of Marconi et a 

D ©04), which are primarily derived 
from bright, high-accretion rate objects. For objects with 
A < rhcri,, we apply the corrections derived from lHol {1999b) 
and Ho et al. (2000), whose samples are essentially all below 
this accretion rate. Although the exact point at which we de- 
termine the divide in bolometric corrections is somewhat ar- 
bitrary, this choice is consistent with theory and observations 
in both cases, and it is not a significant uncertainty as the two 
sets of corrections give similar results for the objects in the 
transition region. 

Observations of X-ray binaries have revealed a number of 
distinct spectral states, such as the high soft state, low hard 
state, quiescent state, etc. (see McClintock & Remillard 2005 
for a review). Moreover, the same black hole can exhibit dif- 
ferent states at different times. A primary cause of the vari- 
ous states is changes in the mode of accretion, as suggested 
by Narayan (1996) and Esin, McClintock & Narayan (1997). 
According to this model, accretion flows with m greater than 
a critical value m cr i t are radiatively efficient and are well- 
described by a thin disk model. These systems correspond to 
the high soft spectral state. However, for m < m cr ; t , the accre- 
tion flow switches to a two-zone state in which the thin disk 
is restricted to radii R greater than a transition radius R tmns , 
while the gas inside R tmns accretes via a radiatively inefficient 
accretion flow (RIAF), also known as an advection-dominated 
accretion flow (ADAF; Narayan & Yi 1994, 1995b, Narayan 
2004). This is the low hard state, or (at very low m) the qui- 
escent state. The transition radius R Vms is only slightly larger 
than the radius of the innermost stable circular orbit when rh 
is just below m cr i t , but it increases to larger values as rh de- 
creases further (see Yuan & Narayan 2004). Consequently, 
the radiative efficiency is fairly close to e° r = 0.1 when rh is 
just below J7i cr it, but it decreases rapidly (but smoothly) as rh 
decreases below m cl -i t (see Equation II II below). 

Observations suggest that the different accretion 
states observed in X -ray b inaries exist as well in AGN 
( iNaravan et alJ fl995l fl996t lMeie3l200!l iMaccarone etafl 
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12001 lYuan & Naravanl kCXM). It appears that the transition 
between accretion states occurs at more or less the same 
critical Eddington ratio m = m cr j t , regardless of the mass 
of the accreting black hole. Although the value of m CI i t is 
somewhat uncertain (see the discussion below on hysteresis), 
observations of black hole binaries ( Maccarona 120031) as 
well as theoretical e x tensions of accretion models (e.g. , 
[ Naravan & YH Il995bl lEsin. McClintock. & Naravai] 119971 
Mever. Liu. & Mever-Hofmeister 2000) suggest m cr ; t ~ 0.01. 
Recent observations of both a bimodal distribution of Ed- 
dington ratios at low redshift ( Marchesini et al. 2004) and the 
distribution of objects in the L-M plane (Jester 2005) also 
sugge st that the transition occurs at m ~ 0.01 (see lCao & Xd 
2005). Note that the similarity between m CI n and niu is 
entirely coincidental. 

In accretion models, the radiative efficiency does not 
depend on the absolute value of M or M, but only on 
the dimensionless accretion rate m (e .g.. IChe n et alJl !995t 
lEsin. McClintock. & Naravad 1 19971 lAbramowiczl 120051) 
Based on observations and theoretical models, we model the 
radiative efficiency e r (m) as 

f e° if m > Wait 

This particular choice for the effici ency factor was orig inally 
suggested by early ADAF models (Naravan & Yi 1995b), but 
it is generally consistent with observations and ensures that 
the radiative efficiency is continuous across the critical Ed- 
dington ratio m cr i t . With this model, we can determine the 
accretion rate m for any black hole given the bolometric lu- 
minosity L and black hole mass M. The model of e r (m) may 
actually overestimate the accretion rate for a given luminos- 
ity at low m <C m cr i t . This is because ADAFs tend to suffer 
large mass loss through winds (Narayan & Yi 1994, 1995a; 
Blandford & Begelman 1999). As a result, m decreases with 
decreasing radius, so that the efficiency depends on whether 
one considers the accretion rate at the black hole or at the 
transition radius. In terms of m at the black hole, which is 
the quantity of interest for this paper, e r decreases less steeply 
with decreasing m, say as ~ (rh/m„it) l l 2 \ i.e., the accretion 
rate for a given luminosity would be less than Equation ( II It 
predicts. However, for our analysis we choose the prescrip- 
tion given in Equation M It . which is the conservative choice 
since it errs on the side of overestimating the accretion rate at 
low A. 

What value of m a -i t should we use? It has been observed 
in X-ray binaries that there is a hysteresis phenomenon in the 
transition between radiatively inefficient and radiatively effi- 
cient accretion flows (see, e.g. Maccarone & Coppi 2003 and 
references therein; Barret & Olive 2002). The transition from 
the radiatively inefficient state (ADAF) to the radiatively ef- 
ficient (thin-disk) state often occurs at a fairly large accre- 
tion rate m ~ 0.1 whereas the reverse transition (from radia- 
tively efficient, to inefficient accretion) occurs at a much lower 
m ~ 0.01. This implies that when m increases with time, the 
critical accretion rate may be as high as »i cl -i t ~ 0.1, whereas 
when m decreases with time m a i t ~ 0.01. Since the peak iof 
quasar activity occurred at high redshift and the activity de- 
clined more or less monotonically to the present epoch, one 
expects that the lower value of m CI i t ~ 0.01 is more relevant 
for our analysis. However, the possibility of hysteresis does 
introduce a degree of uncertainty in our model. Therefore, 
we consider in our subsequent analysis both m cr i t = 0.01 and 



fWciit = 0.1, and ultimately find that the choice does not quali- 
tatively change our conclusions. 

Figure [2] shows the distribution of dimensionless accre- 
tion rates m obtained from the observed sample using Equa- 
tion (II It (upper panels). The left panel shows the result for 
fWciit = 0.1, the right for rii alt = 0.01, essentially bracketing the 
range allowed by the hysteresis effect discussed above. In 
each case, we show the distribution for all objects (empty 
histograms) and for just those objects which are morpho- 
logically identified as ellipticals (morphological type index 
T < -2, shaded histograms). We also show the value of 
mh (dashed line), the accretion rate for which the black hole 
growth timescale is equal to the Hubble time. Although the 
result is slightly stronger for the elliptical subsample, the crit- 
ical point is the same in both cases: even after accounting 
for the fact that the black holes are in a radiatively inefficient 
state and thus the actual accretion rate is significantly larger 
than the observed luminosity Eddington ratio, most local ob- 
jects are accreting at Eddington ratios m <C riin; i.e. at very 
low accretion rates at which significant fractional mass gain 
in less than the age of the Universe is not possible. Although 
the choice of m C rit systematically shifts the inferred m at very 
low m <C m cl -it , it has almost no effect in the range of interest, 
m ~ 10~ 3 - 10 _1 . The bolometric corrections we adopt also 
slightly increase the characteristic m over that implied by e.g. 
a constant eg = 17 (in other words, simpler bolometric cor- 
rections only strengthen this point), although the difference 
is negligible for all but the smallest (m < 10~ 4 ) and largest 
(m > 0.1) accretion rates. 

We can also attempt to use the observed sample accre- 
tion rate distribution to (very roughly) infer that of the pop- 
ulation as a whole. We estimate a simple visibi lity volum e 
correction, using the distances to each object in lHd |2002). 
Comparing the maximum distances or distances within which 
> 90% of each (dynamical or reverberation-mapping mass 
measurements) subsample are located, this gives a rough ra- 
tio of the subsample effective volumes. Alternatively, we 
correct the relative number density of sources in each L 
bin to the value implied by th e z ~ low-lumin osity SDSS 
AGN luminosity function of lHao et akl ll2005l) . In either 
case, the estimate is sim ilar. Whether or not we include the 
Marchesini et al. (2004) sample in this estimate also makes 
little difference. In the lower panels of Figure 13 we use this 
correction (tech nically, the luminos ity function correction in- 
cluding the Marches ini et alj j2004) objects) to estimate a rel- 
ative m distribution (squares), with assumed Poisson errors 
(logarithmically scaled). We caution that these distributions 
are rough estimates and should be taken heuristically, as in 
truth the selection effects involved in measuring individual 
black hole masses are much more complicated than, e.g., a 
simple magnitude limit. We also show the observed SP SS 
Eddington ratio distributions estimated in lYuetalJJ2005ll . in 
a more robust statistical sense, but for relatively high m > 
10~ 3 - 10~ 2 . Again, the normalization of these is essentially 
arbitrary, but the trend with m should be robust. We show 
these observations as the filled dots, for systems with velocity 
dispersions a = 70, 80, 90, 100, HOkms" 1 . 

In all cases, the m distributions broadly agree, and suggest 
that the m distribution can be approximated as a power-law 
(or series of power laws) over a wide range in m. The solid 
lines in the lower panels of Figure |2] show power -laws with 
slopes of —1 /1 .23 estimated from I Yu et alJ ll2005l) and -1 /2 
predicted in Hopkins et al. (2005e), which give reasonable 
approximations over a wide range of m and demonstrate the 
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FIG. 2. — Upper panels: The distribution of dimensionless accretion rate (m = M/Me6<i) calculated through eq. II II and binned by Alogfji = 0.25, from the 
sample of Ho 1 2002). Left panels assume m cr j t = 0.1, right panels assume m aA = 0.01. Shaded histograms show the distribution for AGN with elliptical host 
galaxies only. Lower panels: Estimated relative number density of objects as a function of m (black points), corrected using a simple visibility volume correction 
(alternatively, corrected to match the local AGN luminosity functions of Hao et al. 1 2005 )) for the combined samples of Ho 1 2002 ) and Marchesini et al. 1 2004) 
(squares). Dots show the SDSS Eddington ratio distributions determined in Yu et al. 1 2005) for systems with velocity dispersions a = 70— 1 10km s . The dashed 
line in all panels corresponds to the accretion rate below which the black hole growth timescale M/M is larger than the age of the Universe. 



range in slope allowed. For a power-law rii down to a turnover 
at m m \ n , the differential duty cycle takes the form 

df 



dlog;7t 



= (J3\n\0m f i m )m 



-0 



(12) 



Essentially all of the observational uncertainty regarding duty 
cycles, the proper effective volume and number density cor- 
rections for these samples, and the depth at which the m dis- 
tribution turns over (as compared to e.g. selection incomplete- 
ness) is neatly contained in theparameter m m in- The estimates 
from observations in Figure |2J do at least provide a reason- 
able upper limit m m \ n < 10~ 5 , which allows us to calculate 
the quantities of particular interest (e.g. upper limits to the 
mass gain in low-efficiency phases), given that (3 is also rea- 
sonably well-constrained in the range above. Although these 
estimates are by no means a rigorous determination of the m 
distribution, they do provide approximate constraints and up- 
per limits, which imply that our subsequent calculations (e.g. 
the mass-weighted effective m and e r ) cannot be qualitatively 
changed within the allowed uncertainty in Figure |2] 

2.3. Accretion Rates as a Function of Host Morphology 

Using the intrinsic mass accretion rates estimated above, 
accounting for both the dependence of radiative efficiency on 
accretion rate and spectral shape (bolometric corrections) on 
luminosity and accretion rate, we can consider whether there 
is a difference in the distribution of accretion rates in systems 
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FIG. 3. — The black hole growth times cale ? bh =t s/m as a function of host 
galaxy morphology, for the sample from Ho 1 2002 ). Dashed line marks the 
Hubble time. 



of different morphological types. In Figure [5] we show the 
black hole growth timescale ?bh = ts/rh for each object in the 
sample of iHol J2002I) as a function of the host galaxy mor- 
phology (morphological type in dex), where the host galaxy 
Hubble types are determined in d e Vaucouleurs et all d 1 99 ll) 
(RC3). The dashed line in the figure shows the Hubble time, 



Low-Efficiency Black Hole Accretion 



7 



10' 1 



O O n 



o " o 



10 s 10'° 10" 10 ,s 

Growth Timescale t BH [yr] 



FIG. 4. — The black hole growth timescale as a function of black hole 
mass, open circles the entire sample and filled ellipticals. Dashed line shows 
the Hubble time. 



t H w 14Gyr. It is clear that essentially all objects with moder- 
ate to large accretion rates (i.e. t Bii < f H ) are spirals, whereas 
almost every elliptical has a negligible mass accretion rate 
(tsn ^ ?h)- This is also clear in Figure[2] an d is well known 
from a number of other observations (see e.g., Hec kman et al.l 
2004). 

Observational studies of the z = black hole mass function 
estimated from observations of galaxy luminosity functions 
and the distribution of spheroid velocity dispersions have es- 
tablished that the total present black hole mass density is dom- 
inated by black holes in elliptical hosts, with at most a compa- 
rable contribution from black holes in SO hos t galaxies (e.g. 
[ Aller & Richstond 120021 iMarconi et alJl2004l IShankar et al.1 
2004). If, as is clear in Figure|3] the large majority of present 
elliptical galaxies have ?bh 3> this implies that the majority 
of the present black hole mass density cannot have been ac- 
cumulated in the observed low-accretion rate phase, but that 
ellipticals must have experienced a previous phase with high- 
efficiency accretion. 

Furthermore, regardless of whether or not all bright quasar 
activity is the result of galaxy mergers or interactions, it is dif- 
ficult to imagine a major merger which does not excite bright 
quasar activity (so long as a reasonable amount of cold, ro- 
tationally supported gas is present in either of the merging 
galaxies). If elliptical galaxies are formed from the merger of 
gas-rich spiral progenitors, as is expected in hierarchical sce- 
narios of structure formation (e.g. Toomre & Toomre 1972, 
Toomre 1977; for a review, see Barnes & Hernquist 1992), 
then it is difficult to construct a scenario whereby an elliptical 
galaxy forms without some period of rapid, high-/7i accretion, 
which need only occur for > 10 6 yr to dominate an entire Hub- 
ble time worth of accretion at the present (z = 0) rate. Indeed, 
the gas dissipation that would fuel the growth of supermas- 
sive black holes is required if galaxy mergers are to explain 
the high phase space densities (e.g. Hernquist et al. 1993), 
kinematic properties (e.g. Cox et al. 2005), and fundamental 
scaling relations (e.g. Robertson et al. 2005) of ellipticals. 

3. MODELS FOR THE EVOLUTION OF BLACK HOLE MASS 

3.1. If Black Holes Accrete at Constant m 

Could each black hole accrete at roughly a constant m dur- 
ing its entire lifetime? If this is the case, then objects with 
m <C mh = 0.003 have ?bh 3> fff, and do not gain mass in a 



Hubble time, whereas objects with m > 10~ 2 with ?bh *C ?h 
experience many e-foldings. Thus there should be a strong 
change in the black hole mass distribution at ?bh = with 
much higher black hole masses for ?bh < f# ■ Figure [4] shows 
the black hole mass as a function of growth time fBH for each 
object in our sample (ellipticals filled, late-types open), with 
the dashed line marking the Hubble time. Clearly, the pre- 
dicted trend is not evident; if anything, the most massive ob- 
jects in our sample have the lowest accretion rates (highest 
?b h). This trend has bee n studied in greater statistical detail 
bv lHeckman et al.M2004l) . 

3.2. Models Involving Time Evolution ofm 

The conclusion of the previous subsection is that m must 
vary with time and that the mass growth in many systems must 
occur in an earlier phase of high m accretion. We now ask 
whether low efficiency, low accretion rate growth can give a 
significant fractional contribution to the absolute black hole 
mass density or individual black hole masses at z = 0. The 
answer to this depends on the time history of the accretion 
rate, for which we need to adopt some reasonable model. 

Consider a general model in which m = 1 at early times and 
then decays as a power-law in time, peak phase with m—\, 



m=l/(l+(f/f e )"). 



(13) 



Any such model will naturally produce a power-law distribu- 
tion of accretion rates at low m (Equation II 21 ) with (3 = l/rj 
and m m i n = (f# /tQ)~ v (and in the appropriate limits, can also 
represent exponential or step-function decay). The observed 
m distributions are broadly consistent with the above model 
for ij ~ 1.25-2.0 (se e Figure[2j. Simulations of quasar light 
curves (Hopkins et al. 2005 g) suggest 77 w 2, as does connect- 
ing z = accretion rates at t ~ f# to a bright m ~ 1 phase of 
duration ~ 10 7 - 10 8 yr implied by a wide range of observa- 
tions (see Martini 2004, and references therein). We therefore 
adopt r\ = 2 (this choice also allows analytical solutions be- 
low), but our calculations are only changed by factors ~ rj/2 
for other values of 77. Simulations and theoretical m odels 
(e.g.lSilk & Reesll998[lFabianll999HCiotti & OstrikerKOOlt 
Wvith e & Loebll2003t Pi Matteo et al.ll2003l 120041 120051) . a 
wide range of observations (e.g., Martini 2004), and our esti- 
mated upper limit m m i n < 10~ 5 (simply from the existence of a 
significant population of low-ra observed ellipticals) imply a 
bright (high-?7i) quasar lifetime ~ 10 6 - 10 s yr, setting a fairly 
strict upper limit t q < ts ~ 4 x 10 7 yr. 

The above model of m(t) with r] = 2 gives, for an initial 
black hole mass Mo, the following result for the mass growth 
as a function of time, 

M(f)=M exp|^ tan -1 (j-^j j; (14) 

i.e. the object grows exponentially as M = Moexp(f /ts) for 
t < t Q and asymptotes to M f = Mo<^ (7^/ ' 2,( ' e//fs, as t — > 00. If 
(q ^> ts, then growth is large, but this is because the black 
hole spends a great deal of time in bright, high-;« phases be- 
fore the accretion rate slowly drops off. In fact, even after a 
Hubble time, the accretion rate would still be large. To reach 
m < 0.01 by f = f// requires fg < 3/5, yielding a small mass 
gain at low m; the maximum fraction of the present black 
hole mass which can be accreted at m < 0.01 by the present is 
~8%withr G «1.8f 5 . 

More complex and host-galaxy dependent quasar light 
curves have been obtained in hydrodynamical simulations of 
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galaxy mergers incorporating gas cooling, a multiphase in- 
terstellar medium, Bondi-Hoyle accretion determined from 
the surrounding gas properties, and thermal feedback from 
accreti on |Springel & Hernquisil 120021 120031 ISpringel et alJ 
2005b). Despite the fact that these light curves show 
most objects spend significantly more time at luminosities 
and accretion rates well be low the peak quasar luminos- 
ity JHopkins et al.ll20"05cUel) . the authors find that the mass 
growth is dominated by efficient phases, with > 70 - 80% 
of mass accumulated in bright, optically observable rh ~ 1 
quasar phases (rising to ~ 100% of the mass in the most mas- 
sive objects M > 10 9 M Q ), and the remaining mass essentially 
entirely acquired in t he short following ph ase of declining rh 
at moderate rh ~ 0. 1 ( Hopkins et al l2005q) . 

3.3. The Mass-Weighted Eddington Ratio and Radiative 
Efficiency 

Using the parameterization of the quasar light curve given 
in Equation H3\ . we can define an "effective" mass-weighted 
accretion rate (rh), i.e. 



(rh) = 



1 



m(t)Mdt, 



(15) 



which can be evaluated numerically for a given tg/ts- This 
gives (rh) w 0.50, 0.34 for t g = 10 7 yr and tg = t$ respectively, 
in good agreement with other estimates of the typical accre- 
tion rates at wh ich black holes gain most of their mass (e.g., 
iKollmeier et aIll2005T) . It requires a very large tg > 10 ?s to 
give low (m) < 0.02, but such long quasar lifetimes are ruled 
out by a number of arguments as disc ussed above. Using 
the fitted power laws of lYu et all (12005) gives a significantly 
shallower decay, 77 w 1.23, but similar (rh) w 0.21, 0.14 for 
tg = 10 7 yr and tg = ts- If we use the binned rh distribution of 
Figure|2]directly, we obtain (rh) w 0.3, 0.2 for tg = 10 7 yr and 
tg = ts, with a negligible difference if we assume m cl -i t = 0.01 
or m cr it = 0.1. The difference due to how we estimate the cor- 
rection to an "intrinsic" rh distribution is also small compared 
to the expected Poisson noise from the small sample. Assum- 
ing the distribution cuts off at the lowest observed rh sets an 
upper limit on the duty cycle (lower limit on (m)) independent 
of tg (see Equationll2>. giving (rh) > 0.18. 

We can similarly define an "effective" mass or luminosity- 
weighted radiative efficiency, replacing rh with e r in Equa- 
tion (1151 above. Our "canonical" radiative efficiency e° r = 0.1 
can be factored out of this equation, and we can then numer- 
ically calculate the resulting (e r ) for a given tg/ts and m cl -i t . 
For a relatively large tg = ts, we obtain (e r ) = 0.81 e° (m cr i t = 
0.1); 0.93 e° r (m crit = 0.01). For a shorter timescale t Q ~ 10 7 yr 
these estimates increase to 0.95 e° r , 0.98 ej?. Likewise, if we 
adopt just the observed distribution of accretion rates as the 
PDF for accretion rate over each AGN lifetime as for (rh), we 
obtain (e r ) = 0.92 e° (m clit = 0.1); 0.96 e° (m crit = 0.01) for t Q ~ 
10 7 yr,and (e r ) = 0.72 e° r (^ = 0.1); 0.85 e° r (m cdt = 0.01) for 
tQ ~ ts (and values closer to l.Oe^ for tg < ts). If we were 
to weight by luminosity instead of mass, we would be even 
further biased towards (e r ) ~ e° r . 

Thus, for reasonable values of the quasar lifetime, the "ef- 
fective" radiative efficiencies expected are quite similar to 
the canonical values adopted for bright quasars; essentially, 
this is a restatement of our previous derivation that most 
mass/luminosity is accumulated/radiated in radiatively effi- 
cient, high accretion rate phases. 



3.4. Comparison with Expected Quiescent Accretion Rates 

We can also attempt to estimate the accretion rates and 
growth timescales of the elliptical galaxies observed. If we as- 
sume that the gas in these galaxies is virialized and is in spher- 
ical hydrostatic equilibrium and adopt a Hernquist ( 1990) pro- 
file for the static potential (set by stars and dark matter) we can 
determine the central sound speed and density, and estimate 
the accretion rate at the Bondi-Hoyle rate. Using the M — a 
and M-M sr ,h relations m easured by Trem aine et alJ (120021) 



(16) 



and[Marconi & Huntl (12003). respectively, we obtain 
1 M sph 

<BH —Ol— — — fdyn 



/gas Mbh 



1A 12 /0.01\ (M M . a W a \ 

10 yva {iz){-M-){m^)^ 



ft 

where a ~ 1 is a constant which depends on the profile and 
gas equation of state, w 0.6 for the assumptions above and 
7 = 5/3, fdyn is the spheroid dynamical time, / gas is the gas 
mass fraction, and Mm-ct is the expected M -a black hole 
mass. 

Comparison with Figure|4]suggests that Equation J 16b pro- 
vides a good order-of-magnitude estimate for the ellipticals in 
our observed sample (typical M ~ 10 8 M Q ). Of course, there 
is no particular evidence in the observations for the weak de- 
pendence on a given by our model. But this further suggests 
that these objects must have gone through a previous bright 
phase, to heat the gas to virial temperatures and lower the gas 
fraction to ~ 1% typical in observed ellipticals. Moreover, 
note that ?bh <x Mm-o/M; i.e. if a black hole is undermassive, 
it will take longer to grow, meaning it is unlikely that these 
black holes grew to their observed M -a masses via slow, 
quiescent accretion. 

Finally, in essentially all cases where the Bondi ac- 
cretion rates of l o w-m AGN have been dete r mined 
llFabian & Canizaresl H988L iDi Matteo et alJ l2000t 1200 U 
lLoewenstein et alJl2001tlBower et alJl2003HPellegrinil2005|) 
as well as the same in X-ray binaries (e.g. lPerna et al.1 2003). 
it is observed that objects are actually accreting well below the 
Bondi rate. Thus, the fact that our estimate above provides an 
approximate agreement with the observed growth timescale 
suggests that we may actually be significantly overestimating 
the accretion rates of low-m AGN, making our conclusions 
even stronger. This may be related, for example, to the fact 
that in the Milky Way case, an external observer would likely 
associate a significant amount of nearby but (owing to our 
proximity compared to AGN in other galaxies) separately re- 
solved emission with Sgr A*, because the AGN luminosity is 
not large enough to overwhelm these sources (see the discus- 
sion in lHol2002l) . 

4. CONCLUSIONS 

We have determined the distribution of accretion r ates o f 
z = supermassive black holes from the sample oflHol (12002). 
including the effects of changing radiative efficiency at low 
accretion rates in radiatively inefficient accretion flows and 
the dependence of bolometric corrections on both luminosity 
and accretion rate. We find that most local black holes have 
very low accretion rates rh ~ 10 -4 , well below the minimum 
accretion rate mu = 0.003 required to grow significantly over 
the entire age of the Universe. This is especially true in the 
case of supermassive black holes in elliptical galaxies, imply- 
ing that the dominant part of the integrated black hole mass 
density was formed in radiatively efficient high-m phases, as 
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is expected if these galaxies have undergone mergers in the 
past which would fuel quasar activity and eventually heat or 
expel gas in some form of self-regulated growth (Di Matteo et 
al 2005, Springel et al. 2005a). 

We have further discussed models for the accretion rate as 
a function of time or distribution of accretion rates given the 
observed z = distribution, and find in all cases that the ac- 
cumulated mass should be dominated by mass gained in ef- 
ficient, high-m phases. It is only possible to gain significant 
mass in low-m phases in these models with a high degree of 
fine-tuning of the quasar lifetime, giving values inconsistent 
with both theory and observations. Furthermore, scaling with 
redshift is not expected to change our results significantly, and 
it is increasingly difficult for black holes to gain any mass in 
low-m states if a previous, high-efficiency mass gain period is 
suppressed. 

We do not expect that selection effects will significantly 
chang e ou r resul ts. We have compared directl y the sam- 
ple of iHol (2002) and the co mbined samples of iHol (2002) 
and Marchesini et al. (2004), and find identical qualitative 
results. We have considered separately different classes 
of black holes measured by different techniques, namely 
reverberation-mapping and stellar kinematics, and find that 
the two methods yield similar qualitative results despite hav- 
ing essentially opposite biases. However, because the two 
methods have opposite biases, current measurements with ei- 
ther method only minimally overlap at bolometric luminosi- 
ties L ~ 1O 9 L , allowing the possibility of a bias against 
so urces with accretion rate s m ~ 10 3 — 10 _1 . The addition of 
the M archesini et alJ ll2004l) sample spans this gap in luminos- 
ity, while preserving the deficit of objects at these accretion 
rates. 

Based on the z = distribution of accretion rates, the 
integrated mass of the black hole population must have 
been gained in earlier, radiatively efficient, high-m (m > 
0.1) phases. Likewise, the mass growth of most individ- 
ual black holes must be dominated by high-m rapid growth. 
This is in good agreement with a range of inferred ac- 
cretion histories from observations of quasar p opulations 
(e.g. ISoltanllT98l ISalucci et alJH999t lYu & Tremainell2002t 
Marconi et al. 2004; Shankar et al. 2004), and direct observa- 
tions of quasar accretion rates and their evolution with r edshift 
JVestergaardH2004 iHeckman etaDl2Q04t iMcLure & Dunlopl 
120041 lYu et alJ 120051) . It similarly agrees well with the- 
oretical models of the relation between blac k hole mass 
and host galaxy mass or velocity dispersion {S ilk & Rees 
fl998| pabianlll999l: ICiotti & Ostrikerfl200~fl: IWvithe & Loebl 
2003; Di Matte o et alJ 120 05) and s emi-analytical mod eling 
of ga laxy and quasar formation (e.g. Kauffmann & Haehnelt 
2000). Recently, more detailed modeling of quasar light 
curves and the implied rate of formation of sources from the 
observed luminosity functions (Hopkins et al. 2005a-g) also 
predicts the black hole mass function allowing for luminos- 
ity and host galaxy-dependent quasar lifetimes in which mass 
growth is dominated by a short peak accretion phase of high 
m. Unlike all of these previous observational and theoretical 
constraints, however, our present result is model-independent, 
and does not involve large uncertainties in e.g. the evolution 
of quasars and their fueling mechanisms. 

Because the radiative efficiency of black hole accretion is 
expected and observed to decrease at low m, this argument 
is amplified when considering the radiative luminosity of ac- 



creting black holes. That is, the radiative output is even more 
dominated by high-m phases of accretion. This is consistent 
with the cosmic IR or X-ray backgrounds, synthesis mod- 
els of whic h suggest black hole growth dominated by high-m 
phases (e.g.lComastri et alJl995tlGilli et alJl999HUeda et all 
2003; Hopkins et al. 2005e) and relatively short time spent 
at m rZ JCaol 12005). and also with simple arguments 
which suggest radiative efficie ncies equal to or even larger 
than e,. = 0.1 ( lElvis et all2002tlMerlonl2004l) . Thus, models 
in which quasars spend a very large fraction of their time at 
low luminosities or low accretion rates (e.g. Hopkins et al. 
2005a-g) are completely consistent with both the observed 
Z = distribution of accretion rates and the above observa- 
tions which all suggest that mass growth and radiative output 
are dominated by short, high-m phases. 

Finally, we have not ruled out the possibility that some 
small fraction of individual objects may gain significant mass 
in low-m states. Based on the observed m distribution, this 
fraction is small, < 5%, and we have shown that it does not 
contribute significantly to the integrated black hole mass den- 
sity. However, it could have important implications for the en- 
ergetics of these objects and their relations to their host galax- 
ies, as feedback may couple differently to the surrounding 
medium (in, for example, the form of jets instead of integrated 
disk radiation, given the low-accretion rate state; e.g. Narayan 
& Yi 1995a). Likewise, it is not clear whether the physical 
mechanisms responsible for the M-a relation in most sys- 
tems would apply (or whether alternative mechanisms might 
have the same effect) in low-m systems (although see e.g. 
Churazov et al. 2005). 

To summarize, the observed z = distribution of black hole 
accretion rates constrains the total time black holes spend in 
different states and the relative contribution of these states to 
the integrated final black hole mass and luminosity output. 

(1) Supermassive black holes spend most time in a low- 
efficiency, low accretion rate state. The z = Eddington ra- 
tio distribution implies that the bright quasar lifetime fg < ts- 
This means that ~ 1% of a supermassive black hole's lifetime 
is spent in a bright, high-m quasar phase and ~ 99% is spent 
in a dim low-m phase. 

(2) Nevertheless, the growth of black hole mass is dom- 
inated by the short, high-m phase of evolution. The Ed- 
dington ratio distribution implies an effective mass-weighted 
(m) ~ 0.3, in good agreement with estimates from a wide 
range of other observations of quasars. 

(3) Likewise, the integrated luminosity radiated is domi- 
nated by the high-m, high radiative efficiency growth phase. 
The effective luminosity-weighted radiative efficiency im- 
plied by the observations, coupled with simple models of the 
quasar evolution, is (e r ) > 0.8- 0.9 e° r , where e° ~ 0.1 is the 
radiative efficiency at m = 1 . Again, this is consistent with a 
range of constraints from observations. 

We thus confirm current wisdom regarding points (1) and 
(3), and we obtain a robust result in the case of (2). 
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